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Pulse design without rotating wave approximation 
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We design realizable time-dependent semiclassical pulses to invert the population of a two-level 
system faster than adiabatically when the rotating-wave approximation cannot be applied. Different 
approaches, based on the counterdiabatic method or on invariants, may lead to singularities in the 
pulse functions. Ways to avoid or cancel the singularities are put forward when the pulse spans few 
oscillations. For many oscillations an alternative numerical minimization method is proposed and 
demonstrated. 
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I. INTRODUCTION 

Controlling accurately the internal states of quantum 
two-level systems, realized by real or artificial atoms, 
as in crystal defects, quantum dots, or superconducting 
qubits, is a fundamental task in nuclear magnetic reso¬ 
nance metrology or to develop new quantum technolo¬ 
gies [ll-Q- Pulse engineering is the art and science of 
designing realizable control fields to perform specific op¬ 
erations. We assume here that the field is intense enough 
to be treated semiclassically. The task amounts to solv¬ 
ing an “inverse problem”, as the aim is to determine a 
realizable Hamiltonian that drives the system to specific 
states at a given final time, e.g., the ones that would be 
reached in an ideal adiabatic driving M- The design 
and implementation of the pulse are easier for adiabatic 
dynamics, as the final atomic state is quite insensitive to 
smooth deviations from the ideal pulse, but for faster- 
than-adiabatic processes, the design and its implemen¬ 
tation become more demanding. In addition, when the 
rotating-wave approximation (RWA) m fails for strong 
drivings, the inversion process becomes more involved. 
However, this regime is key to achieve fast control of 
two-level systems, e.g., in nitrogen-vacancy centers El- 

Ill. 

In this paper we focus on speeding up “rapid adia¬ 
batic passage” (RAP) ^pulation inversion processes in 
two-level systems E3 EE El , beyond the RWA. With¬ 
out the RWA, a consistency condition between the diag¬ 
onal and the non-diagonal elements of the interaction- 
picture Hamiltonian, which involves the phase ip{t) of 
the field and its derivative with respect to time, 
must be satisfied. We shall first apply the counterdia¬ 
batic (CD) method ElEEH- In this method, a refer¬ 
ence Hamiltonian, i.e., the Hamiltonian for the two-level 
system without the RWA in the interaction picture, H{t), 
is complemented (or even substituted) by a CD Hamil¬ 
tonian, so that the system follows exactly the 

adiabatic dynamics of Hit). In the dynamics driven by 
H{t) A- Hi{t), the CD term suppresses transitions in the 
instantaneous eigenbasis of H, but allows for transitions 
in the instantaneous eigenbasis of H A Hi. This Hamil¬ 
tonian, however, does not necessarily satisfy the consis¬ 


tency condition mentioned above, unless an appropriate 
rearrangement of terms is performed, and a new unitary 
transformation between interaction and Schrodinger pic¬ 
tures, different from the one for H, is implemented [2l|. 
In section H we show that the consistent field implies in 
general singularities in the Rabi frequency, but the field 
itself is not singular. 

We use as well invariants of motion (2^. Designing 
the invariant is equivalent to imposing the desired dy¬ 
namics. This is easy within the RWA by setting the 
time dependence of independent auxiliary variables cor¬ 
responding to polar, 9{t), and azimuthal angles, /3(t), 
that characterize the invariant eigenstates on the Bloch 
sphere EE- From these angles the time dependences of 
the Hamiltonian components, the Rabi frequency VInit) 
and the detuning A(t), and thus the physical fields may 
be deduced E3. By contrast, without applying the RWA, 
a naive independent design leads to singularities in the 
field. In section HI we show two different ways to avoid 
this problem. When the pulse duration spans only a few 
field oscillations the angles 9 and /3 may be set to can¬ 
cel all singularities and produce a smooth, finite-intensity 
pulse. Instead, if many field oscillations occur, a numer¬ 
ical minimization method to find optimal parameters in 
a predetermined pulse form may be used to invert the 
population. 

We shall set first the basic concepts and notation. As¬ 
suming a semiclassical interaction between the electric 
field, E{t) = Ao(t) cos[v3(t)], and the two-level atom, the 
Hamiltonian of the atom in the Schrodinger picture, in 
the electric dipole approximation, is, see e.g. (23l |. 

H,(t) = ^{a;o(t)(|e)(e|-|5)(5l) 

+ D«(t)[(|e)(5| + |5)(e|)(e*‘^(‘) + e-*^(*))]},(l) 

in the bare basis of the atom \g) = ( (, ), |e) = ( “ ), 
where uJo{t)/(2 tt) is the transition frequency of the atom, 
which may in general depend on time, controlled, e.g., 
by Stark shifts; H/j(t) is the Rabi frequency, assumed 
real (without loss of generality for transitions without 
change in the magnetic number, see EH); and ifit) is the 
time dependent phase of the electric field of the pulse, 
where (f{t) / (27r) is the instantaneous field frequency. The 
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“exact” Hamiltonian, i.e., without applying the RWA, in 
a field-adapted interaction picture is given by [2^ 

H{t) = Ul{Hs - H^)U^, (2) 

where 

HAt) = Y^\e){e\-\g){9\) (3) 

and 

U^{t) = in 

= e-*^(‘)/2|e)(e|+e*^(‘)/2|5)(5l (4) 

is the unitary operator of the transformation. Thus, 


method to speed up an adiabatic population inversion 
process for a two-level systems beyond the RWA, where 
the (reference) Hamiltonian of the system is given by Eq. 
m, or in diagonal form as H{t) = J2n \''T‘{i))^n{t){n{t)\, 
where \n{t)) = |±(t)). The inversion is from |4'(0)) = |g) 
to |4'(t/)) = |e), up to a phase factor, where t = 0 and 
t = tf are the initial and final times of the process, and 
|4'(t)) is the general state of the system. We consider 
a constant field (angular) frequency, ip{t) = to l, so that 
ip(t) = uj^t, and a time dependent transition frequency, 
LUolt) /(27r), that will determine the detuning. Then, from 
Eqs. ® and ([7]), 

nit) = Hfl(l + e-2*-^‘), (9) 

A(t) = uJoit)-ujL- (10) 


Hit) 


ft / -A n\ 
2\n* A j ’ 


with 


H(t)=Hji(t)[l + e-2*‘^«], 


and detuning 


A(t) = ujoit) - (fit). 


(5) 

( 6 ) 

(7) 


The CD term is in general given by, see e.g. [^ . 
Hiit) = ihy^i\dtn)in\ - {n\dtn)\n){n\) 

n 

\m)im\dtH\n){n\ 


= E E ■ 

m^n 


En - E„ 


For the exact two-level system, using Eq. dS]), i/i is given 
by HI 


The eigenenergies of H are E±it) = ±ft£o(t)/2, with 
eoit) = \J A(t)2 -|- |H(t)p, and the (time dependent) 
eigenstates, |±) = |±(t)), are 


l±) 


-(ATgp) 

n* 


\g) + \e) 


1 



(Areoli’ 


( 8 ) 


which satisfy 77(t)|±(t)) = i?±(t)|±(<)). 

The exact two-level system Hamiltonian in the inter¬ 
action picture entails ^pit) in the detuning, i.e., in the 
diagonal elements, see Eq. and its integral, </5(t), in 
the non-diagonal elements, see Eq. ®. By “consistency 
condition” we mean that the elements of a physically al¬ 
lowed interaction picture Hamiltonian must comply with 
the structure set in Eqs. and ©. In particular, 
this structure must be satisfied when designing pulses 
to carry out faster-than-adiabatic inversion processes, as 
we shall see in sections H and HI. Suppose for example 
that the functions H/{(t) and uJoit) are given. Then, not 
any Hamiltonian is allowed (consistent) since the diago¬ 
nal and non-diagonal parts must depend on ipit) and its 
derivative consistently. 


II. THE COUNTERDIABATIC METHOD 

The counterdiabatic approach adds a counterdiabatic 
or “CD” term to some reference Hamiltonian to make 
the exact dynamics adiabatic with respect to the ref¬ 
erence Hamiltonian [T| [TtI - I^ . The formal construc¬ 
tion of the CD term is explained in the original refer¬ 
ences [ISElBil. We shall apply the counterdiabatic 


Hiit) 


ih f -Bi/2 Ai \ 
2Ci -Tit Bi/2 ) ’ 


where Ai(t) = DA — AO, i?i(f) = D*D — tin* is purely 
imaginary, and C'i(t) = A^-|-|Dp. Thus, the total Hamil¬ 
tonian provided by the CD method is 


H + Hi 


ft / -A n\ 
2\ n* Ay’ 


( 11 ) 


where 


n = n + iAi/Ci, (12) 

A = A + iBi/2Ci. (13) 


From Eq. m, the detuning of the total Hamiltonian, 
given by Eq. (US, becomes 


A = iuJo + iBi/2Ci)-uJL- (14) 


The term u}o+iBi/2Ci may be interpreted as a new time- 
dependent transition (angular) frequency, and Wi/(27r) 
is, as in the reference Hamiltonian, the constant field 
frequency. Then, the non-diagonal element D, given by 
Eq. ([12]), should be expressed as D = D_r( 1 -I- 
compare with Eq. ([HI) , in order to satisfy the consistency 
condition, where D/j = n^it) is a new Rabi frequency 
corresponding to the pulse associated with H + Hi. 

To see whether H + Hi is indeed consistent we first 
rewrite the complex D in Eq. m in a convenient form, 

n = Dfl(l-t e“2*‘7) (15) 

= 2Dfl cos(ye“*‘7, (16) 
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where flu and if = if(t) are real, if is given by 


— arg(fi) + 27rn, (17) 


where n is an integer chosen to make f continuous. Once 
f is determined, taking into account Eq. flu is cal¬ 
culated from Eqs. (HID or (HI. However, the consistency 
condition is generally not satisfied, i.e., f ^ f = ujLt. For 
the Allen-Eberly protocol [HI, 


A{t) 


Om sinh 

2SHo 


Tr{t-tfl2) 


tanh 


2to 
Tr{t-tf/2) 

2tn 


(18) 


and parameters Om = 27r x 3 MHz, 6 = 2Tr x 200 MHz, 
to = 0.05 ns, wl = 27r X 10 GHz, and tf = 0.4 ns. Fig. 
[1] (a) shows that f{t) = and if{t) do not coincide. 
These parameters are chosen so that H{t) does not in¬ 
vert the populations of the bare basis, Pg{t) = K5|4>(f))p 
and Pe{t) = |(e|4'(f))p. Fig. [T](b) shows that the Hamil¬ 
tonian H + Hi does invert the population. The CD term 
works formally, but the Hamiltonian H+Hi does not cor¬ 
respond to a field with frequency ljJl/{2tt). If we apply 
the same transformation that relates the Schrodinger pic¬ 
ture Hamiltonian and the interaction-picture Hamil¬ 
tonian H, 


UUt) = e-“-‘/2|e)(e| + (19) 


for f = Wit, following Eq. the Schrodinger picture 
Hamiltonian corresponding to H + Hi does not take the 
form of Eq. o with modified functions for the transition 
and Rabi frequencies. This procedure is schematized in 
Fig. [2]by the boxes around S, /, and S", where S rep¬ 
resents the initial Schrodinger picture driven by Hg and 
I the interaction picture driven by H with given by 
Eq. (HU, the unitary operator that generates the trans¬ 
formation. I' represents the interaction picture with the 
addition of an extra term, H+Hi , and S' the correspond¬ 
ing Schrodinger picture mediated again by Eq. m- 

We might as well assume that Hi is a Hamiltonian 
corresponding to an independent (second) pulse comple¬ 
menting the pulse for H. This, however, leads to the 
same result as interpreting H + Hi as a single pulse, 
since H is negligible versus Hi, i.e., |H| <C \iAi/Ci\ and 
A ^ iBij{2Ci), for the given parameters. The consis¬ 
tency condition is again not satisfied. 

To enforce the consistency condition between the di¬ 
agonal and non-diagonal elements in Eq. (El) we may 
rewrite A, see Eq. d, with the required structure, i.e., 
A = Cjo{t) — if, where wo(t)/(27r) is a new time-dependent 
transition frequency of the atom. Equating this expres¬ 
sion to Eq. d gives 


Ao(l) — Wo — u!l + iBi/2Ci + If, 


which is depicted in Fig. [I](c) with uJo{t) obtained from 
Eq. (fTUl) and A given by Eq. (fTSl) . The unitary transfor¬ 
mation Uf, where f is substituted by if in Eqs. m and 





FIG. 1: (Color online) (a) if{t) = lolI (green solid line) and 
fit) given by Eq. (11711 (blue short-dashed line), (b) Pg driven 
by H (red dashed line) and by H+Hi (bine solid line), and (c) 
ioo{t) = A + If divided by 27r, for the Allen-Eberly protocol 
and parameters: Dm = 27r x 3 MHz, 6 = 2n x 200 MHz, 
to — 0.05 ns, ljl = 27r x 10 GHz, and tf = 0.4 ns. 


o, leads to a different Schrodinger picture S" [2l[, see 
Fig. O The Hamiltonian in the Schrodinger picture S" 
is 


Hs"{t) = 


—Wo 2Hi{ cos f 
2flji cos If Wo 


This procedure solves the consistency condition, but the 
new Rabi frequency, given by Eqs. m or (HU), is singu¬ 
lar at the zeros of cos ip. However, the field, proportional 
to cos If, is well behaved, see Fig. |3l 

In the following section we work out a different ap¬ 
proach based on invariants of the motion. 
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FIG. 2: (Color online) Schematic relation between differ¬ 
ent Schrodinger and interaction picture dynamical equations. 
Each node may represent the dynamical equations, and also 
the Hamiltonian or the states. The rectangular boxes enclose 
nodes that represent the same underlying physics. The solid 
lines represent unitary relations for the linked states and the 
dashed line represents a non-unitary addition of a term to the 
Hamiltonian. 
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FIG. 3: cos (p/(27r) for parameters as in the caption of 

Fig.m 


III. INVARIANT-BASED INVERSE 
ENGINEERING 


Associated with H (t) there are hermitian dynamical 
invariants, I{t), satisfying the invariance condition [^ . 


dl 1 


that may be parametrized as 


I{t) 


filo f COS 9 sin0e \ 

2 \ sin — cos 9 J ’ 


( 20 ) 


( 21 ) 


where Iq is an arbitrary constant angular frequency to 
keep I{t) with dimensions of energy, and 9 = 9{t) and 
P = P{t) are the polar and azimuthal angles in the Bloch 
sphere, respectively. The eigenvalue equation for I(t) 
is I[t)\(f)^^{t)) = A±|(/>4(t)), where the eigenvalues are 


\± = and, consistently with orthogonality and 

normalization, we can choose the basis of eigenstates 




( e -iS- 

/ COS 2 

1 ' 9 

\ sin|e 2 


, I'/'-W) 


• B -iS- \ 

sin^e 2 \ 

9 ■ 

— cos|e 2 J 


( 22 ) 

The solution of the time-dependent Schrodinger equa¬ 
tion, ihdt\'^{t)) = iJ(t)|'F(t)), can be expressed as 

14/(0) = c+e*^+«|0i(0) +C-e*^-W|0'_(0), 


up to a global phase factor, where the c± are time- 
independent coefficients and the 7±(0 are the Lewis- 
Riesenfeld phases. 


7±(i) = - Hit') (l)i{t')'jdt'. 


The Lewis-Riesenfeld phase becomes a global phase if the 
dynamics is carried out by one eigenstate of the invariant 
only. This will be the case in the inversions discussed here 
and leads to an important simplification: 7 ±(t) can be 
ignored to engineer the Hamiltonian. 

From the invariance condition o, with I given by 
Eq. (1^ and H given by Eq. (O, we find 


9 = —2V, n cos (psm{l3 — ip), 

/3 = — A — 2H/{Cot0cos(^cos(/3 — (^). (23) 

The first equation bounds 9 between —2Vii and 20i{. 
These equations are much simplified when the RWA may 
be applied 0- The two-level system Hamiltonian under 
the RWA is, see for example [lOl| . 

HRWAit) = J ( ) 3 


where pit) is absent in the non-diagonal elements of 
HRWAit): compare with Hit) in Eq. ([5]). Now the spher¬ 
ical angles satisfy 0, [H, 0 


rwa[1 


-Hfl;sin/3, 

—A — Vr cot 9 cos p. 


(24) 


The “direct problem” is to solve the systems of differ¬ 
ential equations (l23l) or (l24l) for 9 and /? when Vr and 
A are given, once we fix wq or p for the system (l23ll . see 
also Eq. ©• Instead, in the “inverse problem” we have 
in principle to construct the functions Vr and A from 
9 and /3. When the RWA holds, from the system (l24l) . 
the inversion reduces to simple expressions for Vr and A 
0 , 


RWA I 


—9/ sin/3, 

—P + 9 cot 9 cot p. 


Without the RWA, from the system (051) . taking into 
account Eq. 0, we get 


d 


9 

2 cos p sin a ’ 

—Wo + ^cot 9 cot a. 


(25) 

(26) 
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FIG. 4: (Color online) Interpolated functions (a) 6 and (b) a 
for a population inversion process with parameters wz, = 27r x 
500 MHz and t/ = 5 ns. (c) and (d) show the corresponding 
flH/(27r), given by Eq. (1^ . and A/(27r), given by Eq. (1^ . 
(e) shows the populations of the bare basis, Pg (blue long- 
dashed line) and Pe (red solid line), for these Hamiltonian 
functions, and (f) shows the time-dependent atom frequency, 
uJo/(2n), derived from A in Eq. dTj as cJo(t) = A(t) +ujl- 


where a{t) = j3{t) — ip{t) is introduced to simplify the ex¬ 
pressions. If the inversion strategy is to consider u;o(i) as 
given, which in particular could be constant, and 0(t) is 
designed, the differential equation for a(t), see Eq. (|M)) . 
is problematic, as an arbitrary choice of 9{t) and wo(t) 
will typically lead to singularities on the right hand side, 
and a{i) will, in general, introduce singularities in VLfi{t). 
Eq. (1251) makes clear that a finite, smooth fl/j and a 
phase Lp increasing along many field cycles, require many 
zeros of 0 to compensate for zeros in the denominator. A 
naive choice for 9 is thus doomed to fail. A different ap¬ 
proach is needed. It is useful to rewrite Eq. (051) . taking 
into account Eq. o , as 


carry the dynamics, see Eq. (I22L with 9{t) going from 
0 to TT. In addition, H{t) and I{t) should commute at 
t = 0 and t = tf so that the initial and final states, 
l^^(O)) and are eigenstates of the initial and 

final Hamiltonians. This implies vanishing derivatives of 
9{t) at the initial and final times. We have in summary 
the boundary conditions 

0(O)=0(O) = 0(t/)=O, 9{tf)=^. (28) 


A. Pulse with few field oscillations 

For pulses containing few field oscillations and assum¬ 
ing again a constant (angular) frequency for the pulse 
field, ip = ujl, and ip(t) = uiLt, we may construct 9 so 
that 9 cancels the zeros of cos(/7 in Eq. (I25L and then 
a to compensate the singularities of cot 9 in Eq. (1221), 
bounded by 0 < a < tt not to introduce new singular¬ 
ities in Eqs. (1251) and (l27)) . In the example of Fig. |4] 
we use wz, = Stt X 500 MHz and tf = 5 ns, and interpo¬ 
late 9{t) and a{t) with polynomials 9 = 
a = J2ti=o 9 is constructed to satisfy the boundary 
conditions (l28ll . and to make 9 zero at the five intermedi¬ 
ate zeros of cos ip. 0(1 ns) = 2, 0(1.6 ns) = 2.4, 0(2.5ns) = 
2.8, 0(4 ns) = 2.8, and 0(4.5 ns) = 3 are also imposed to 
force a smooth ascent of 0 up to tt, see Fig. |4](a). At the 
boundary times, tf, = 0,t/, the conditions given by Eq. 
(l28l) imply that limt-,.^,, (0 cot 0 cot a) = —2d(tb), see Eq. 

The polynomial ansatz for a is constructed so that 
a becomes 7r/2 at the singularities of cot0, in this case 
at t = 0 and t = tf. At the boundary times we choose 
as an example d(0) = d(t/) = 0, which corresponds to 
wo(0) = uJo{tf) = 0. The function is additionally tamed 
for smoothness by imposing a{tf/2) = 2, see Fig. |4](b). 
Hr and A calculated from Eqs. (1251) and (l27ll are shown 
in Figs. 111(c) and (d), respectively. The populations of 
the bare basis, Pg{t) and Pe{t), are shown in Fig. |4](e). 
Fig. l4](f) shows the time-dependent transition frequency 
of the atom, uJo{t)/{2Tr), given by Eq. (lH), which changes 
sign. This is possible by varyin g la ser intensities around 
a “light-induced level crossing” [27 1. 


B. Pulse with many field oscillations 


A = — (yj-b d)-I-0cot 0cot Of. (27) 

Even if 0 and a avoid the singularities for a given p in Eq. 
dm), the required consistency between Eqs. d?]) and dm 
(they must be equal) will in general fail if uJo{t) is also 
given. Thus, the proposed strategy is to fix first p{t), 
then design 9{t) to avoid the zeros of cosp in Eq. (ESI) , 
and from there design a(t) to compensate singularities of 
cot 0 in Eq. dm- This produces a smooth Hr given by 
Eq. dm, A follows from Eq. (l27l) . Finally, Wo(t) is 
deduced consistently from Eq. d3- 

For a population inversion from |4'(0)) = |g) to 
l^(^/)) = I®), to a phase factor, we may use \(j)+{t)) to 


For many field oscillations in the pulse the singularities 
to avoid in equations (ESI) and (im become too numer¬ 
ous to apply the previous approach. A simple inversion 
method is to assume sensible specific forms with free pa¬ 
rameters for the functions A and Hr. An example is a 
population inversion process carried out by a linear de¬ 
tuning and a Gaussian Rabi frequency, 

A(t) = a{t-tf/2), (29) 

HR{t) = Hoexp[-A{t-tf/2)^], 

with two free parameters, a and Hq. For a constant tran¬ 
sition (angular) frequency, ujq, taking into account Eqs. 
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FIG. 5: (Color online) Pg driven by HfiwA{t) (blue dashed 
line) and by the exact Hamiltonian H(t) (red solid line) for (a) 
the reference parameters a = {2n)^ x 254.648 MHz^ and flo = 
27r X 2 GHz, and for (b) the optimized parameters a = (Stt)^ x 
272.824 MHz^ and flo = 2tv x 2.202 GHz, for a population 
inversion process. In both cases tf — 0.1 /rs, oiq = 27r x 5 
GHz, and A = ( 271 )^ x 506.606 MHz^. 

© and (l29t . setting (/^(O) = 7r/2 and integrating, the 
phase of the pulse is given by 

ip(t) = —ai^ 12 + (wo + atf /2)t + -k 12. 

Setting /3(0) = 0, we may solve for 9{t) and /3(t) in the 
system (1231) . and minimize — tt]^ numerically to 

determine a and flo- This is a simple alternative to a 
more sophisticated optimal-control-theory approach [ 2 ^ 
or bang-bang methods (^ . 

In a numerical example we first set the reference 
parameters, tf = 0.1 ^s, wq = 27r x 5 GHz, A = 
(27r)2 X 506.606 MHz^, a = {2tt)^ x 254.648 MHz^, and 
Ho = 27r x 2 GHz, for which the Hamiltonian within 
and without the RWA give similar dynamics with unsuc¬ 
cessful population inversions, see Fig. [S](a). For these 
values of tf, ujq, and A, and from these seed parameters 
a and Hq a minimization algorithm provides optimized 
parameters for the population inversion process with the 
exact Hamiltonian i?(t), a = (27r)^ x 272.824 MHz^ and 
Hq = 27r x 2.202 GHz. The same optimized parameters 


do not invert the population when the RWA is applied, 
see Fig. [5] (b). 


IV. DISCUSSION AND CONCLUSIONS 

For a two-level system in a classical field the diago¬ 
nal and non-diagonal elements in the exact interaction- 
picture Hamiltonian H{t) must depend consistently on 
the phase of the field, ^{t), and its derivative, ip{t). This 
makes the inverse engineering methods more complicated 
than with the Hamiltonian within the rotating-wave ap¬ 
proximation, HfiwA{t). Simple attempts using invari¬ 
ants or counterdiabatic methods to implement faster- 
than-adiabatic processes may not satisfy the consistency 
condition in H(t) or lead to singularities. Different ways 
have been shown to circumvent the difficulties. While we 
have presented simple proof-of-principle examples, the 
methods may be adapted to minimize effects of noise and 
decoherence due to the flexibility of the inversion [30l| . 

An open problem is to extend the approaches to sys¬ 
tems where more levels have to be considered. In many 
systems the failure of the RWA is associated with the 
need to include further levels in the theoretical treatment. 
In trapped ions, for example, when the vibrational RWA 
is not applied, and vibrational counter-rotating terms are 
taken into account, the energy levels are distorted and 
the sideband resonances are shifted [3l|. This may be 
understood as a vibrational Bloch-Siegert effect or, equiv¬ 
alently, as the result of Stark shifts of the levels due to 
off-resonant transitions [3l|. Nevertheless, it is possible 
to describe the subspace of the two states in an isolated 
anticrossing by an approximate 2x2 Hamiltonian that 
takes into account the effect of further levels perturba- 
tively by means of a level shift operator [^. The 
current approaches could then be applied but the details 
are left for a separate study. 
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